Tidal deformability of neutron stars with realistic equations of state and their 
gravitational wave signatures in binary inspiral 
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The early part of the gravitational wave signal of binary neutron star inspirals can potentially yield 
robust information on the nuclear equation of state. The influence of a star's internal structure on 
the waveform is characterized by a single parameter: the tidal deformability A, which measures the 
star's quadrupole deformation in response to the companion's perturbing tidal field. We calculate 
A for a wide range of equations of state and find that the value of A spans an order of magnitude 
for the range of equation of state models considered. 

An analysis of the feasibility of discriminating between neutron star equations of state with 
gravitational wave observations of the early part of the inspiral reveals that the measurement error 
in A increases steeply with the total mass of the binary. Comparing the errors with the expected 
range of A, we find that Advanced LIGO observations of binaries at a distance of 100 Mpc will probe 
only unusually stiff equations of state, while the proposed Einstein Telescope is likely to see a clean 
tidal signature. 

PACS numbers: 04.40. Dg, 26.60.Kp, 97.60. Jd, 95.85. Sz 



I. INTRODUCTION AND SUMMARY 

The observation of inspiraling binary neutron stars 
(NSs) with ground-based gravitational-wave detectors 
such as LIGO and Virgo may provide significantly more 
information about neutron-star structure, and the highly 
uncertain equation of state (EOS) of neutron-star matter, 
than is currently available. The available electromagnetic 
observations of neutron stars provide weak constraints 
from properties such as the star's mass, spin, and grav- 
itational redshift (see for example [TJ [3] ) . Simultaneous 
measurements of both the mass and radius of a neutron 
star [31 HI |5J [7] , on the other hand, have the potential 
to make significantly stronger constraints. These mea- 
surements, however, depend on detailed modeling of the 
radiation mechanisms at the neutron-star surface and ab- 
sorption in the interstellar medium, and are subject to 
systematic uncertainties. 

Another possibility for obtaining information about 
the neutron star EOS is from the inspiral of binary neu- 
tron stars due to gravitational radiation. The tidal dis- 
tortion of neutron stars in a binary system links the EOS 
describing neutron-star matter to the gravitational-wave 
emission during the inspiral. Initial estimates showed 
that for LIGO, tidal effects change the phase evolution 
only at the end of inspiral, and that point particle post- 
Newtonian waveforms can be used for template-based 
detection [5J O HO]. With the projected sensitivities of 
later-generation detectors, however, effects which can be 
neglected for the purpose of detection may become mea- 
surable in the strongest observed signals. 

While EOS effects are largest during the late inspiral 
and merger of two neutron stars where numerical sim- 



ulations must be used to predict the signal, Flanagan 
and Hinderer showed that a small but clean tidal sig- 
nature arises in the inspiral below 400 Hz This 
signature amounts to a phase correction which can be 
described in terms of a single EOS-dependcnt tidal de- 
formability parameter A, namely the ratio of each star's 
induced quadrupole to the tidal field of its companion. 
The parameter A depends on the EOS via both the NS 
radius R and a dimensionless quantity k 2 , called the Love 
number [121 [T3J []J] : A = 2/(3G)k 2 R 5 . 

The relativistic Love numbers of polytropic 1 EOS were 
examined first by Flanagan and Hinderer [TTJ [TS] and 
later by others in more detail [THl EZ]. Flanagan and 
Hinderer also examined the measurability of the tidal 
deformability of polytropes and suggested that Advanced 
LIGO could start to place interesting constraints on A 
for nearby events. However, they used incorrect values 
for &2) which overestimated A by a factor of ~ 2 — 3 
and were therefore overly optimistic about the potential 
measurability. In addition, polytropes are known to be 
a poor approximation to the neutron star equation of 
state, and there may be significant differences in the tidal 



Polytropes are often written in two forms. The first form is 
expressed as p = Ke 1+1 / n , where p is the pressure, e is the 
energy density, K is a pressure constant, and n is the polytropic 
index. The second form, is given by p = Kp 1+1 / n , where p 
is the rest-mass density, defined as the baryon number density 
times the baryon rest mass. The first form was mainly used in 
the recent papers 1151 1161 [T7l . However, the second form is more 
commonly used in the neutron-star literature and is more closely 
tied to the thermodynamics of a Fermi gas. We will use both 
forms as was done in Ref. |16| . 
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deformability between polytropes and "realistic" EOS. 
In this paper, we calculate the deformability for realistic 
EOS, and show that a tidal signature is actually only 
marginally detectable with Advanced LIGO. 

In Sec. [II] we describe how the Love number and tidal 
deformability can be calculated for tabulated EOS. We 
use the equations for k 2 developed in [15j . which arise 
from a linear perturbation of the Oppenheimer-Volkoff 
(OV) equations of hydrostatic equilibrium. In Sec. Ill we 
then calculate k 2 and A as a function of mass for several 
EOS commonly found in the literature. We find that, 
in contrast to the Love number, the tidal deformability 
has a wide range of values, spanning roughly an order of 
magnitude over the observed mass range of neutron stars 
in binary systems. 

As discussed above, the direct practical importance of 
the stars' tidal deformability for gravitational wave ob- 
servations of NS binary inspirals is that it encodes the 
EOS influence on the waveform's phase evolution during 
the early portion of the signal, where it is accurately mod- 
eled by post-Newtonian (PN) methods. In this regime, 
the influence of tidal effects is only a small correction to 
the point-mass dynamics. However, when the signal is 
integrated against a theoretical waveform template over 
many cycles, even a small contribution to the phase evo- 
lution can be detected and could give information about 
the NS structure. 



Following [llj, we calculate in Sec. IV the measurabil- 
ity of the tidal deformability for a wide range of equal- 
and unequal- mass binaries, covering the entire expected 
range of NS masses and EOS, and with proposed detector 
sensitivity curves for second- and third- generation detec- 
tors. We show that the measurability of A is quite sensi- 
tive to the total mass of the system, with very low-mass 
neutron stars contributing significant phase corrections 
that are optimistically detectable in Advanced LIGO, 
while larger-mass neutron stars are more difficult to dis- 
tinguish from the k 2 = case of black holes [TSJ [T7] . In 
third-generation detectors, however, the tenfold increase 
in sensitivity allows a finer discrimination between equa- 
tions of state leading to potential measurability of a large 
portion of proposed EOSs over most of the expected neu- 
tron star mass range. 

We conclude by briefly considering how the errors 
could be improved with a more careful analysis of the 
detectors and extension of the understanding of EOS ef- 
fects to higher frequencies. 

Finally, in the Appendix we compute the leading or- 
der EOS-dependent corrections to our model of the tidal 
effect and derive explicit expressions for the resulting cor- 
rections to the waveform's phase evolution, extending the 
analysis of Ref . [11] . Estimates of the size of the phase 
corrections show that the main source of error are post- 
1 Newtonian corrections to the Newtonian tidal effect 
itself, which are approximately twice as large as other, 
EOS-dependent corrections at a frequency of 450 Hz. 
Since these point-particle corrections do not depend on 
unknown NS physics, they can easily be incorporated into 



the analysis. A derivation of the explicit post-Newtonian 
correction terms is the subject of Ref. [18] . 
Conventions: We set G = c = 1. 



II. CALCULATION OF THE LOVE NUMBER 
AND TIDAL DEFORMABILITY 



As in [TT] and [15], we consider a static, spherically 
symmetric star, placed in a static external quadrupolar 
tidal field £ij. To linear order, we define the tidal de- 
formability A relating the star's induced quadrupole mo- 
ment Qij to the external tidal field, 



Qij A^j . 

The coefficient A is related to the I 
tidal Love number k 2 by 
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= 2 dimcnsionless 
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The star's quadrupole moment Qij and the external 
tidal field £ij are defined to be coefficients in an asymp- 
totic expansion of the total metric at large distances r 
from the star. This expansion includes, for the met- 
ric component g u in asymptotically Cartesian, mass- 
centered coordinates, the standard gravitational poten- 
tial m/r, plus two leading order terms arising from the 
perturbation, one describing an external tidal field grow- 
ing with r 2 and one describing the resulting tidal distor- 
tion decaying with r~ 3 : 



(1 + fltt) 



3 Qij 



(3) 



where n % — x % jr and Qij and £ij are both symmetric and 
traceless. The relative size of these multipolc components 
of the perturbed spacetime gives the constant A relating 
the quadrupole deformation to the external tidal field as 
in Eq. 0. 

To compute the metric ([3]), we use the method dis- 
cussed in |15j . We consider the problem of a linear static 
perturbation expanded in spherical harmonics following 
[TT?] . Without loss of generality we can set the azimuthal 
number m = 0, as the tidal deformation will be axisym- 
metric around the line connecting the two stars which 
we take as the axis for the spherical harmonic decompo- 
sition. Since we will be interested in applications to the 
early stage of binary inspirals, we will also specialize to 
the leading order for tidal effects, I = 2. 

Introducing a linear I — 2 perturbation onto the spher- 
ically symmetric star results in a static (zero- frequency) , 
even-parity perturbation of the metric, which in the 
Regge- Wheeler gauge [20] can be simplified [15] to give 



ds 1 



_ e 2*W [l + H(r)Y 20 (9,^]dt 2 
+e 2A « [l-H{r)Y 20 {6, V )]dr 2 
+r 2 [1 - K(r)Y 20 (9, <p)] [dO 2 + sin 2 9d^ 2 ) 



(4) 
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where K(r) is related to H{r) by K'{r) = H'(r) + 
2H(r)<&'(r). Here primes denote derivatives with re- 
spect to r. The corresponding perturbations of the per- 
fect fluid stress-energy tensor components are 5T ° = 
-6e(r)Y 20 (6, tp) and 5Tj = 5p(r)Y 20 {6, <p), where e is the 
energy density and p the pressure. The function H{r) 
satisfies the differential equation 



fip 2A 3 
.^__ 2 ($') 2 + 2$" + ^A' 

- 7 $'-2$'A'+ Lftf + A')) H 
r r I 



$' - A' ) H' + H" = 0. 



(5) 



Here / is given by 



5e = fSp 



(6) 



which for slow changes in matter configurations corre- 
sponds to / = de/ dp. 

The method of calculating the tidal perturbation for a 
general equation of state table is similar to the method of 
calculating moment of inertia in the slow rotation approx- 
imation 21j. The specific implementation we use follows 
the moment of inertia calculation in Appendix A of |2J , 
via an augmentation of the OV system of equations 2 : 
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The second-order differential equation for H is separated 
into a first-order system of ODEs in terms of the usual 
OV quantities m r 3 , p(r), and e(p), as well as the addi- 
tional functions H(r), (3{r) = dH/dr, and the equation 



of state function f(p) (recall / = de/dp): 
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These are combined with Eqs. ([7|-pTj|, and the aug- 
mented system is solved simultaneously. The system is 
integrated outward starting just outside the center us- 
ing the expansions H{r) = a$r 2 and /3(r) = 2aor as 
r — * 0. The constant a® determines how much the 
star is deformed and can be chosen arbitrarily as it can- 
cels in the expression for the Love number. The ODE 
for H(r) outside the star, where T^„ = 0, has a gen- 
eral solution in terms of associated Legendre functions 
Q\ (r/m — 1) ~ r~ 3 at large r, and Pf ( r / m — 1) ^ r 2 
at large r. The boundary conditions that determine the 
unique choice of this solution follow from matching the 
interior and exterior solutions and their first derivatives 
at the boundary of the star, where r = R. By compari- 
son with Eq. ([3| , the coefficients of the external solution 
can then be identified with the axisymmetric tidal field 
and quadrupole moment via £Y2o(0,(p) 



QY 20 (6,<p) = Qi 



Ei^rtn 3 , and 



-\£i 



as was done in |15j . 



Here, £ and Q are the magnitudes of the 1 = 2, m = 
spherical harmonic coefficients of the tidal tensor and 
quadrupole moment respectively. 
Defining the quantity 



RP(R) 
H(R) 



(13) 



for the internal solution, the / = 2 Love number is 

k 2 = ^(l-2C) 2 [2 + 2C(y-l)-y] 

<2C[6-3y + 3C(5y-8)} 

4C 3 [13 - lly + C(3y - 2) + 2C* 2 (1 + y)} 
3(1 - 2C) 2 [2 -y + 2C(y - 1)] ln(l - 2C) 



(14) 



2 Here we present the equations in terms of the radial coordinate 
r; the extension to the enthalpy variable r) used in [2] is straight- 
forward. 

3 We use m r for the mass enclosed within radius r instead of m(r) 
to avoid confusion with the total mass of the star, which we will 
label m. 



where C = mj R is the compactness of the star. 

For stars with a nonzero density at the surface (for 
example strange quark matter or an incompressible n = 
polytrope), the term (//r)($' + A') in Eq. ^ blows up 
at the surface r = R and H'(r) is no longer continuous 
across the surface. Following the discussion in [55] for 
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an n — polytrope, this discontinuity leads to an extra 
term in the expression above for y: 

_ R(3(R) 47ri? 3 e_ 
V ~ H(R) Tn 

where e_ is the density just inside the surface. 



(15) 



III. LOVE NUMBERS AND TIDAL 
DEFORMABILITIES FOR CANDIDATE EOS 

Differences between candidate EOS can have a signif- 
icant effect on the tidal interactions of neutron stars. In 
this paper we consider a sample of EOS from Refs. [21125] 
with a variety of generation methods and particle species. 
The sample is chosen to include EOS with the largest 
range of behaviors for k 2 {m/R\ k<i(rn) and A(m) rather 
than to fairly represent the different generation methods. 
We also restrict ourselves to stars with a maximum mass 
greater than 1.5 M©, which is conservatively low given 
recent neutron-star mass observations [2H 1251 1251 [271 128j . 
We consider 7 EOS with just normal npep matter 
(SLY [29 , API and AP3 [30], FPS [31], MPA1 [32], 
MSI and MS2 [33]), 8 EOS that also incorporate some 
combination of hyperons, pion condensates, and quarks 
(PS [34], BGN1H1 [35], GNH3 [36], HI and H4 [37], 
PCL2 [35], ALF1 and ALF2 [35]), and 3 self-bound 
strange quark matter EOS (SQM1-3 [3D]). A brief de- 
scription of these EOS and their properties can be found 
in [2] [23]. 

The generic behavior of the Love number k 2 is shown in 
the top panel of Fig. [I] as a function of compactness m/R 
for different types of EOS. The two types of polytropes, 
energy and rest-mass density polytropes, are shown in 
gray. They coincide in the limit m/R — > where e — ► p 
as the star's density goes to zero, and in the limit n — > 
where e(p) and p{p) are both constant. This can be seen 
from the first law of thermodynamics, 



d— = —pd-, 
P P 



(16) 



which relates e to p. 

The sequences labeled "Normal" correspond to the 15 
EOS with a standard nuclear matter crust, and the 3 se- 
quences labeled "SQM" correspond to the crustless EOS 
SQM1-3 where the pressure is zero below a few times 
nuclear density. Within these two classes, there is little 
variation in behavior, so we do not explicitly label each 
candidate EOS. 

The bottom panel of Fig. [I] shows k 2 (m) for the realis- 
tic EOS, which is more astrophysically relevant because 
mass, not compactness, is the measurable quantity dur- 
ing binary inspiral. Unlike the quantity k2(jn/R), k2{m) 
depends on the constant K for polytropes, so polytropic 
EOS are not shown. There is more variation in k 2 for 
fixed mass than for fixed compactness. 

The behavior of these curves can be understood as fol- 
lows: The Love number k 2 measures how easily the bulk 
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FIG. 1: Top panel: Love number as a function of com- 
pactness. Gray dotted curves are energy density polytropes 
(p = A"e 1+1 / n ), and gray solid curves are rest-mass den- 
sity polytropes (p = Kp 1+1 ^"). Both polytropes are the 
same for n = 0. EOS with only npep matter are solid and 
those that also incorporate 7r/hyperon/quark matter are dot- 
dashed. The three SQM EOS are dashed and overlap. They 
approach the n — curve at low compactness, where hi has 
a maximum value of 0.75 as m/R — > 0. Bottom panel: Love 
number as a function of mass for the same set of realistic EOS. 
Note that there is more variation in hi between different EOS 
for fixed mass than for fixed compactness. 



of the matter in a star is deformed. If most of the star's 
mass is concentrated at the center (centrally condensed) , 
the tidal deformation will be smaller. For polytropes, 
matter with a higher polytropic index n is softer and 
more compressible, so these polytropes are more centrally 
condensed. As a result, k 2 decreases as n increases. The 
limiting case n — represents a uniform density star and 
has the largest Love number possible. The Love num- 
ber also decreases with increasing compactness, and from 
Eq. (14 1 it can be seen that k 2 vanishes at the compact- 
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ness of a black hole (m/R= 0.5) regardless of the EOS 
dependent quantity y [16j [17] . 

Normal matter EOS behave approximately as poly- 
tropes for large compactness. However, for smaller com- 
pactness, the softer crust becomes a greater fraction of 
the star, so the star is more centrally condensed and &2 
smaller. For strange quark matter, the EOS is extremely 
stiff near the minimum density, and the star behaves ap- 
proximately as an n — polytrope for small compact- 
ness. As the central density and compactness increase, 
the softer part of the EOS has a larger effect, and the 
star becomes more centrally condensed. 

The parameter that is directly measurable by gravi- 
tational wave observations of a binary neutron star in- 
spiral is proportional to the tidal deformability A, which 
is shown for each candidate EOS in Fig. [2] The values 
of A for the candidate EOS show a much wider range of 
behaviors than for &2 because A is proportional to fc2-R 5 , 
and the candidate EOS produce a wide range of radii 
(9.4-15.5 km for a 1.4 M star for normal EOS and 8.9- 
10.9 km for the SQM EOS). See Table [I] 



TABLE I: Properties of a 1.4 M Q neutron star for the 18 EOS 
discussed in the text. 



EOS 


R(km) 


m/R 


fc 2 


A(10 36 gcm 2 s 2 ) 


SLY 


11.74 


0.176 0.0763 


1.70 


API 


9.36 


0.221 


0.0512 


0.368 


AP3 


12.09 


0.171 


0.0858 


2.22 


FPS 


10.85 


0.191 


0.0663 


1.00 


MPA1 


12.47 


0.166 


0.0924 


2.79 


MSI 


14.92 


0.139 


0.110 


8.15 


MS2 


13.71 


0.151 


0.0883 


4.28 


PS 


15.47 


0.134 


0.104 


9.19 


BGN1H1 


12.90 


0.160 


0.0868 


3.10 


GNH3 


14.20 


0.146 


0.0867 


5.01 


HI 


12.86 


0.161 


0.0738 


2.59 


HI 


13.76 


0.150 


0.104 


5.13 


PCL2 


11.76 


0.176 


0.0577 


1.30 


ALF1 


9.90 


0.209 


0.0541 


0.513 


ALF2 


13.19 


0.157 


0.107 


4.28 


SQM1 


8.86 


0.233 


0.098 


0.536 


SQM2 


10.03 


0.206 


0.136 


1.38 


SQM3 


10.87 


0.190 


0.166 


2.52 



For normal matter, A becomes large for stars near the 
minimum mass configuration at roughly 0.1 Mq because 
they have a large radius. For masses in the expected 
mass range for binary inspirals, there are several differ- 
ences between EOS with only npe/j, matter and those 
with condensates. EOS with condensates have, on aver- 
age, a larger A, primarily because they have, on average, 
larger radii. The quark hybrid EOS ALF1 with a small 
radius (9.9 km for a 1.4 M Q star) and the nuclear matter 
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FIG. 2: Tidal deformability A of a single neutron star as a 
function of neutron-star mass for a range of realistic EOS. The 
top figure shows EOS that only include npefj, matter; the mid- 
dle figure shows EOS that also incorporate 7r/hyperon/quark 
matter; the bottom figure shows strange quark matter EOS. 
The dashed lines between the various shaded regions repre- 
sent the expected uncertainties in measuring A for an equal- 
mass binary inspiral at a distance of D = 100 Mpc as it passes 
through the gravitational wave frequency range 10 Hz-450 Hz. 
Observations with Advanced LIGO will be sensitive to A in 
the unshaded region, while the Einstein Telescope will be able 
to measure A in the unshaded and light shaded regions. See 
text below. 
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only EOSs MSI and MS2 with large radii (14.9 km and 
14.5 km, respectively, at 1.4 M & ) are exceptions to this 
trend. 

For strange quark matter stars, there is no minimum 
mass, so the radius (and therefore A) approaches zero as 
the mass approaches zero. At larger masses, the tidal 
deformability of SQM stars remains smaller than most 
normal matter stars because, despite having large Love 
numbers, the radii of SQM stars are typically smaller. 

Error estimates A A for an equal-mass binary inspiral 
at 100 Mpc are also shown in Fig. [2] for both Advanced 
LIGO and the Einstein Telescope. They will be discussed 
in the next section. 



IV. MEASURING EFFECTS ON 
GRAVITATIONAL RADIATION 
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We wish to calculate the contribution from realistic 
tidal effects to the phase evolution and resulting grav- 
itational wave spectrum of an inspiraling neutron star 
binary. In the secular limit, where the orbital period is 
much shorter than the gravitational radiation reaction 
timescale, we consider the tidal contribution to the en- 
ergy E and energy flux dE/dt for a quasi-circular inspi- 
ral using the formalism developed by Flanagan and Hin- 
derer which adds the following leading-order terms 
to the post-Newtonian point-particle corrections (PN-PP 
corr.): 



E(x) 



E(x) 



1 



-Mrjx 



1 + (PN-PP corr.) 



9^4x 5 + 1 ~2 



32 



-rj 2 x 5 



1 + (PN-PP corr.) 



mi M 5 



(17) 



(18) 



Here Ai = A(mi) and A2 = A(m 2 ) are the tidal deforma- 
bilities of stars 1 and 2, respectively. M = mi + 7772 is 
the total mass, 7/ = toitt^/M 2 is the dimensionless re- 
duced mass, and x is the post-Newtonian dimensionless 
parameter given by x = (cjM) 2 / 3 , where ui is the orbital 
angular frequency. One can then use 



dx/dt = 



E 



dE/da 



(19) 



to estimate the evolution of the quadrupole gravitational 
wave phase <I> via d$/dt = 2u — 2x 3 / 2 /M. 

Each equation of state gives in this approximation a 
known phase contribution as a function of mi and m2, 
or as a function of the total mass M — mi + 771,2 and 
the mass ratio 7772/77*1, via A(mi) and A(ttt,2) for that 
EOS. Although we calculated A for individual neutron 
stars, the universality of the neutron star core equation 
of state allows us to predict the tidal phase contribution 



FIG. 3: Weighted A for a range of chirp mass M and di- 
mensionless reduced mass 77, for three of the EOSs considered 
above. The values of 77 equal to {0.25, 0.242, 0.222} cor- 
respond to the mass ratios 7712/7711 = {1.0, 0.7, 0.5}. Also 
plotted (as in Fig. [5| are the uncertainties AA in measuring 
A for a binary at 100 Mpc between 10 Hz-450 Hz. The solid, 
dashed, and dotted curves correspond to AA for 77 = 0.25, 
0.242, and 0.222 respectively. 



for a given binary system from each EOS. Following [IT] , 
we discuss the constraint on the weighted average 



A = 



2G 



777l + 127772 w n m 2 

-A(mi) + 



12m 



mi 



777 2 



-A(m 2 ) 



, (20) 



which reduces to A in the equal mass case. The contri- 
bution to d^/dx from the tidal deformation, which adds 
linearly to the known PP phase evolution, is 



d$ 

dx 



195 r^A 
8 M 5 r] ' 



(21) 



The weighted average A is plotted as a function of chirp 
mass M = (mi Jn2 ) 3/5 /M 1/5 in Fig. |] for three of the 
EOS and for three values of 77: equal mass (77 = 0.25), 
large but plausible mass ratio (41] (77 = 0.242), and ex- 
tremely large mass ratio (77 = 0.222). 

We can determine the significance of the tidal effect 
on gravitational waveforms in a given frequency range 
by considering the resulting change in phase accumulated 
as a function of frequency. In the case of template-based 
searches, for example, a drift in phase of half a cycle leads 
to destructive interference between the signal and tem- 
plate, halting the accumulation of signal to noise ratio. 
The phase contributions to binary neutron stars of vari- 
ous masses from a range of realistic tidal deformabilities 
are plotted in Fig. |4j 

The post-Newtonian formalism itself is sensitive to 
high-order corrections at the frequencies at which the 
tidal effect becomes significant; as reference, we show in 
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FIG. 4: The reduction in accumulated gravitational wave 
phase due to tidal effects, $ 3 .5,pp(/gw) - $3.5,a(/gw), is 
plotted with thick lines as a function of gravitational wave 
frequency, for a range of A appropriate for realistic neutron 
star EOS and the masses considered. The 3.5 post-Newtonian 
TaylorT4 PN specification is used as the point-particle refer- 
ence for the phase calculations. For reference, the difference 
in accumulated phase between 3.0 and 3.5 post-Newtonian 
orders of T4 (thin dashed line), and the difference between 
3.5 post-Newtonian T4 and 3.5 post-Newtonian Tl (thin dot- 
ted line) are also shown. Phase accumulations are integrated 
from a starting frequency of 10 Hz. 



expansions, as well as that from varying the form of the 



post-Newtonian Taylor expansion from T4 to Tl. 4 An 
accurate knowledge of the underlying point-particle dy- 
namics will be important to resolve the effects of tidal 
deformation on the gravitational wave phase evolution 
at these frequencies. 

The half-cycle or more contribution to the gravita- 
tional wave phase at relatively low frequencies suggests 
that this effect could be measurable. Flanagan and Hin- 
derer |llj first calculated the measurability for frequen- 
cies below 400 Hz, where the approximations leading 
to the tidal phase correction are well-justified. We ex- 
tend the same computation of measurability to a range 
of masses and mass ratios. We take noise curves from 
the projected NS-NS optimized Advanced LIGO config- 
uration |33], as well as a proposed noise spectrum of the 
Einstein Telescope [45 . These noise curves are repre- 
sentative of the anticipated sensitivities of the two detec- 
tors. Our results do not change significantly for alternate 
configurations which have similar sensitivities in the fre- 
quency range of interest. 

We also extend the computation to a slightly higher 
cutoff frequency. As estimated in the Appendix, our cal- 
culation should still be fairly robust at 450 Hz, as the 
contributions to the phase evolution from various higher 
order effects are O(10%) of the leading order tidal contri- 
bution. The uncertainty in the phase contribution from 
a given EOS is therefore significantly smaller than the 
order of magnitude range of phase contributions over the 
full set of realistic EOS. 

The rms uncertainty AA in the measurement of A is 
computed using the standard Fisher matrix formalism 
[4"5] . Assuming a strong signal h and Gaussian detector 
noise, the signal parameters l have probability distribu- 
tion p (9 l ) cx exp (-(1/2) TijS&SQi), where 56* = 9 l - § l 
is the difference between the parameters and their best- 
fit values 6 l and Fy = (dh/dO* , dh/88 1 ) is the Fisher 
information matrix. The parentheses denote the inner 
product defined in [46J. The rms measurement error in 
9 l is given by a diagonal element of the inverse Fisher, or 
covariance, matrix: AO" 1 — J (r -1 )". 

Using the stationary phase approximation and neglect- 
ing post-Newtonian corrections to the amplitude, the 
Fourier transform of the waveform for spinning point 
masses is given by h(f) = Af~ 7 ^ 6 exp(W), where the 
point-mass contribution to the phase 'F is given to 3.5 
post-Newtonian order in Ref. |47j . The tidal term 



<5* 



tidal 



117Ax 5 / 2 
8r?M 5 



(22) 



4 For an explanation of the differences between T4 and Tl, see [421 

Fig.|4]the phase difference between the 3.0PN and 3.5PN |23]. 
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obtained from Eq. (A. 5) adds linearly to this, yielding 
a phase model with 7 parameters (t c , <f> c , Ai, 77, f3, a, A), 
where (3 and a are spin parameters. We incorporate 
the maximum spin constraint for the NSs by assuming 
a Gaussian prior for (3 and a as in |46j . The uncertain- 
ties computed will depend on the choice of point-particle 
phase evolution, but we assume this to be exactly the 
3.5PN form for the current analysis. 

The rms measurement uncertainty of A, along with the 
uncertainties in chirp mass M. and dimensionless reduced 
mass 77, are given in Table [TT| and plotted in Figs. [2] and [3] 
from a single-detector observation of a binary at 100 Mpc 
distance with amplitude averaged over inclinations and 
sky positions. If the best-fit A is zero, this represents a 
l-<7 upper bound on the physical A. A signal with best- 
fit A > A A would allow a measurement rather than a 
constraint of A, with 1-a uncertainty of AA. 

We obtain the following approximate formula for the 
rms measurement uncertainty AA, which is accurate 
to better than 4% for the range of masses 0.1 M < 
mi, m2 < 3.0 Mq and cutoff frequencies 400 Hz < / cn d < 
500 Hz: 



AA 



M 



m 2 
mi 



fend 

LIz~ 



D 



lOOMpc 



(23) 

where a — 1 .0 x 10 42 g cm 2 s 2 for a single Advanced LIGO 
detector and a = 8.4 x 10 40 gem 2 s 2 for a single Einstein 
Telescope detector. 

Our results show that the measurability of tidal effects 
decreases steeply with the total mass of the binary. Esti- 
mates of the measurement uncertainty for an equal-mass 
binary inspiral in a single detector with projected sen- 
sitivities of Advanced LIGO and the Einstein Telescope, 
at a volume-averaged distance of 100 Mpc and using only 
the portion of the signal between 10 — 450 Hz, are shown 
in Fig. [2] together with the values of A predicted by var- 
ious EOS models. Measurability is less sensitive to mass 
ratio, as seen in Fig. [3j Comparing the magnitude of the 
resulting upper bounds on A with the expected range for 
realistic EOS, we find that the predicted A are greatest 
and the measurement uncertainty AA is smallest for neu- 
tron stars at the low end of the expected mass range for 
NS-NS inspirals of (1 M Q - 1.7 M ) [48]. 

In a single Advanced LIGO detector, only extremely 
stiff EOS could be constrained with a typical 100 Mpc 
observation. However, a rare nearby event could allow 
more interesting constraints, as the uncertainty scales as 
the distance to the source. Rate estimates for detec- 
tion of binary neutron stars are often given in terms of a 
minimum signal-to-noise p c = 8; a recent estimate [49] is 
between 2 and 64 binary neutron star detections per year 
for a single Advanced LIGO interferometer with a volume 
averaged range of 187 Mpc. The rate of binaries with a 
volume averaged distance smaller than 100 Mpc trans- 
lates to roue hly (100/187) J ~ 15% of this total detection 
rate, but over multiple years of observation a rare event 
could give measurements of A with uncertainties smaller 



than the values in Table |H| (e.g. with half the tabled un- 
certainty at 1.9% the total NS-NS rate). 

Using information from a network of N detectors with 
the same sensitivity decreases the measurement uncer- 
tainty by approximately a factor of 1/y/N |50j . giving 
more reason for optimism. However, we should also note 
that, in some ways, our estimates of uncertainty are al- 
ready too optimistic. First, AA only represents a 68% 
confidence in the measurement; a 2AA error bar would 
give a more reasonable 95% confidence. In addition, our 
Fisher matrix estimates are likely to somewhat underesti- 
mate the measurement uncertainty in real non-Gaussian 
noise. 

In contrast to Advanced LIGO, an Einstein Telescope 
detector with currently projected noise would be sensitive 
to tidal effects for typical binaries, using only the signal 
below 450 Hz at 100 Mpc. The tidal signal in this regime 
would provide a clean signature of the neutron star core 
equation of state. However, an accurate understanding 
of the underlying point-particle phase evolution is still 
important to confidently distinguish EOS effects. 

TABLE II: The rms measurement error in various binary pa- 
rameters (chirp mass M, dimensionless reduced mass 77, and 
weighted average A of the tidal deformabilities) for a range of 
total mass M and mass ratio mj/mi, together with the signal 
to noise ratio p, using only the information in the portion of 
the inspiral signal between 10 Hz < / < 450 Hz. The dis- 
tance is set at 100 Mpc, and the amplitude is averaged over 
sky position and relative inclination. 

Advanced LIGO 



M (Mq) 


m 2 /m 1 AM/M A77/77 AA(10 


36 gem 2 s 2 ) 


P 


2.0 


1.0 


0.00028 0.073 


8.4 


27 


2.8 


1.0 


0.00037 0.055 


19.3 


35 


3.4 


1.0 


0.00046 0.047 


31.3 


41 


2.0 


0.7 


0.00026 0.058 


8.2 


26 


2.8 


0.7 


0.00027 0.058 


18.9 


35 


3.4 


0.7 


0.00028 0.055 


30.5 


41 


2.8 


0.5 


0.00037 0.06 


17.8 


33 



Einstein Telescope 



M (M Q ) m 2 /m 1 AM/M Ari/rj AA(10 36 gem 2 s 2 ) p 



2.0 
2.8 
3.4 
2.0 
2.8 
3.4 
2.8 



1.0 
1.0 
1.0 
0.7 
0.7 
0.7 
0.5 



0.000015 0.0058 
0.000021 0.0043 
0.000025 0.0038 
0.000015 0.0058 
0.000021 0.0045 
0.000025 0.0038 
0.000020 0.0048 



0.70 
1.60 
2.58 
0.68 
1.56 
2.52 
1.46 



354 
469 
552 
349 
462 
543 
442 



Expected measurement uncertainty will decrease if we 
can extend the calculation later into the inspiral. From 
Eq. (23), A A at 500 Hz is approximately 79% of its value 



at 450 Hz. The dominant source of error in the tidal phas- 
ing at these frequencies are post-Newtonian effects which 



9 



scale as Xx 7 / 2 and do not depend on any additional EOS 
parameters. These terms are computed in Ref. |18j . and 
when they are incorporated into the analysis, the result- 
ing phase evolution model can be used at slightly higher 
frequencies. These terms also add - 10%(//450 Hz) 2 / 3 
to the strength of the tidal signature. 

Higher-order tidal effects and nonlinear hydrodynamic 
couplings, which depend on unknown NS microphysics, 
are smaller than post-Newtonian effects by factors of ~ x 
and ~ x 2 , so they become important later in the inspi- 
ral, where the adiabatic approximation that the mode 
frequency is large compared to the orbital frequency also 
breaks down. At this point we can no longer measure 
only A, but an EOS- dependent combination of effects 
including higher multipolcs, nonlincarity, and tidal reso- 
nances. 

However, information in the late inspiral could also 
constrain the underlying neutron-star EOS. Read et 
al. [51] estimated potential measurability of EOS effects 
in the last few orbits of binary inspiral, where the grav- 
itational wave frequency is above 500 Hz, using full nu- 
merical simulations. The EOS used for the simulation 
was systematically varied by shifting the pressure in the 
core while keeping the crust fixed. The resulting models 
were parameterized, either by a fiducial pressure or by 
the radius of the isolated NS model, and measurability 
in Advanced LIGO was estimated. Such numerical simu- 
lations include all the higher order EOS effects described 
above, but the I = 2 tidal deformability parameter A 
should remain the dominant source of EOS-dependent 
modification of the phase evolution. We therefore expect 
it to be a better choice for a single parameter to charac- 
terize EOS effects on the late inspiral. 

The numerically simulated models of |51j can be re- 
parameterized by the A of the 1.35 Mq neutron stars 
considered 5 . The uncertainty of measurement for the 
new parameter A can be estimated from Tables II- V of 
[5T] , In the broadband Advanced LIGO configuration of 
Table IV, it is between 0.3 and 4xl0 36 gcm 2 s 2 for an 
optimally oriented 100 Mpc binary, or between 0.7 and 
9xl0 36 gcm 2 s 2 averaged over sky position and orienta- 
tion. However, in the NS-NS optimized LIGO config- 
uration of Table III, which is most similar to the Ad- 
vanced LIGO configuration considered in this paper, the 
expected measurement uncertainty is more than several 
times A for all models. These estimates should be consid- 
ered order-of-magnitude, as numerical simulation errors 
are significant, and the discrete sampling of a param- 
eter space allows only a coarse measurability estimate 
which neglects parameter correlations. In contrast to 
the perturbative/post-Newtonian estimate of EOS effects 
calculated in this paper, EOS information in the signal 



before the start of numerical simulations is neglected. 
The estimate is complementary to the measurability be- 
low 450 Hz estimated in this paper. 



V. CONCLUSION 

We have calculated the relativistic / = 2 Love number 
&2 and resulting tidal deformability A for a wide range 
of realistic EOS in addition to polytropes. These EOS 
have tidal deformabilities that differ by up to an order 
of magnitude in the mass range relevant for binary neu- 
tron stars. However, the estimated uncertainty AA for 
a binary neutron star inspiral at 100 Mpc using the Ad- 
vanced LIGO sensitivity below 450 Hz is greater than 
the largest values of A except for very low-mass binaries. 
The uncertainty for the Einstein Telescope, on the other 
hand, is approximately an order of magnitude smaller 
than for Advanced LIGO, and a measurement of A will 
rule out a significant fraction of the EOS. 

Advanced LIGO can place a constraint on the space of 
possible EOS by obtaining a 95% confidence upper limit 

can also 



IV 



The piecewise polytrope EOS {2H, H, HB, B, 2B} have Ai. 35Ms 
of {0.588, 1.343, 1.964, 2.828, 10.842} X 10 36 gem 2 s 2 , respec- 
tively. 



of \{M,rj) < 2AX(M,rj). The tables in Sec._ 
be scaled as follows: For a network of N detectors the 
uncertainty scales roughly as AX/y/N, and for a closer 
signal we have AA(£>/100 Mpc). 
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APPENDIX: ACCURACY OF THE PHASING 
MODEL 



To assess the accuracy of the simple phase evolution 
model, we compute the corrections to the tidal phase 
perturbation due to several EOS-dependent effects: the 
leading order finite mode-frequency terms, higher order 
tidal effects, and nonlinear hydrodynamic couplings. For 
simplicity, we will only derive the phase corrections for 
one star with internal degrees of freedom coupled to a 
point mass. The terms for the other star simply add. 
For such a binary system, the Lagrangian can then be 
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written as 



L 



1 , 1 , , nM 2 
-jjMr 2 + -T]Mr 2 ip 2 + 

2 2 r 



1 



1 



1 1 / . 2 

~7.Qijk^ijk i ri\ 2~ ( QijkQijk ^^QijkQijk 



^3 QijQ jkQki' 



(A.l) 



deformability constant A3 is denned by Qijk — —^^ijk- 
The quantities luq and CJ03 are the ^ = 2 and I = 3 /-mode 
frequencies, and a is a coupling constant for the leading 
order nonlinear hydrodynamic interactions. In general, 
one would need to sum over the contributions from all 
the modes, but other modes contribute negligibly in the 
regime of interest for the above model (see [55]). Post- 
Newtonian effects on the Lagrangian for the binary are 
derived in Ref. [TH] and can simply be added to those 
derived here. 



Here, the star's static mass quadrupole Qij parameter- 
izes the I — 2 modes of the star, which can be treated as 
harmonic oscillators that are driven below their resonant 
frequency by the companion's tidal field. The tensor Qijk 
parameterizes the star's mass octupole degrees of free- 
dom, and £ij and Sijk are the / = 2 and I = 3 tidal tensors 
respectively, which are given by £ij = didj{— m 2 /r) and 
£ijk — didjdk(~Tn2/r) in Newtonian gravity. The I = 3 



We will be interested in finding an effective descrip- 
tion of the dynamics of the system for quasi-circular in- 
spirals in the adiabatic limit, where the radiation reac- 
tion timescale is long compared to the orbital timescale. 
From equilibrium solutions to the Euler-Lagrange equa- 
tions derived from this Lagrangian, the following radius- 
frequency relation is obtained: 



1 



3Am 2 w 10 / 3 9Ato 2 w 10 / 3 uj 2 



20A 3 m 2 w 14 / 3 18am^ 16 / 3 
M 7 /3 TOl M 8 / 3 m 1 



27A 2 m§w 2 °/3 
M 10 / 3 m 2 



(A.2) 



M 5 /3 mi M 5 / 3 mi uj 2 

The equilibrium energy, obtained by reversing the signs of the potential energy terms in the Lagrangian, is given by: 



E = -- V M 5 /W/ 3 
2 ' 



1 



9Am 2 w 10 / 3 45Ato 2 w 10 / 3 uj 2 65A 3 to 2 w 14 / 3 



M 5 / 3 



mi 



M 5 / 3 r 



M 7 /3 mi 



60am|w 16 / 3 
M 8 / 3 mi 



63A 2 m^ 2 °/ 3 
M 10 / 3 m 2 



(A.3) 



The correction to the energy flux E = — \{QijQij): where Qfj — iLr 2 {ri Lr n? — \&ij) + Qij is the total quadrupole 
moment, is 



E = 



32 



^10/3^10/3 



6Aw 10 / 3 f n m 2 
M 2 / 3 



mi 



M 



12Aw 10 / 3 uj 2 



M 2 / 3 



mi u)q 



80A 3 m 2 w 14 / 3 



M 7 / 3 



mi 



36am 2 w 16 / 3 /_to 2 



9A 2 w 2 °/ 3 



M 5 / 3 mi V 2 M +1 y ' M 4 / 3 m? \" U l/~ 



2m 2 6toj 



(A.4) 



Using the formula cP^/duj 2 — 2(dE/dui)/E in the stationary phase approximation and integrating twice leads to the 
final expression for the tidal phase correction: 



<5# = - 



9Ax 5 / 2 (m l + 12m 2 \ 45Ax 5 / 2 uj 2 /8mi + 155to 2 \ 125 A32: 9 / 2 m 2 



IQrjM 5 \ mi 
135am 2 a; 11 / 2 / mi 
352t]M 8 \ ro7 



UO877M 5 w 2 



mi 



12 ?7M 7 mi 



18m 2 



3A 2 a; 5 ( M 2 - 2m 2 M - 83m 2 . 



64r/M 10 V 
I 



(A.5) 



We will analyze the information contained in the por- 
tion of the signal at frequencies / < 450 Hz. This is 
slightly higher than previously considered, and we now 
argue that in this frequency band, the simple model of 
the phase correction is still sufficiently accurate for our 
purposes. We will evaluate all of the corrections for the 
case of equal masses m\ — m 2 = m. An estimate of 
the fractional errors for the case of to = 1.4 M© and 



R = 15 km is given in parentheses. 

1. Post-l-Newtonian corrections (~ 10%). 

These corrections give rise to terms oc Ax 7 / 2 that 
add to those in Eq. ( A.5 1 . The explicit form of these 
terms is computed in Ref. [18] and they depend on 
the NS physics only via the same parameter A as 
the Newtonian tidal terms, so they can easily be in- 
corporated into the data analysis method. Prclimi- 
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nary estimates indicate that for equal masses, these 
post-1 Newtonian effects will increase the tidal sig- 
nal. 

Adiabatic approximations ( < 1%). 
The approximation that the radiation reaction time 
is much longer than the orbital time is extremely 
accurate, to better than 1%; see Fig. 2 of Ref. [TT] . 
which compares the phase error obtained from 
numerically integrating the equations of motion 
supplemented with the leading order gravitational 
wave dissipation terms to that obtained analyti- 
cally using the adiabatic approximation. 
The accuracy of the approximation w<Cwo can be 



estimated from the fractional correction to (A.5l, 
which is - (815/1144)(w/w ) 2 ~ 0.18(//7oT, 
where f — uj/tt and /o = ujq/(2t:). For typical 
NS models the I — 2 /-mode frequency is [S3] 



kHz 



0.78 + 1.64 



1.4 M, 



^3 



1/2 



R 



10 km 



-3/2 



(A.6) 



so that the fractional correction is ~ 0.012 for 
/ = 450 Hz and for a conservatively low /-mode 
frequency of /o = 1700 Hz. 

Higher order tidal effects (< 0.7%). 
The I = 3 correction to the gravitational wave 
phase (A. 5 1 is smaller than the / = 2 contribu- 
tion by a factor of - (25 / '351) (k 3 /k 2 )(m/R)- 2 x 2 ~ 
0.007, for m/R = 0.14 and a stiff n = 0.5 poly- 
trope. Here, we have defined the I = 3 Love number 
fc 3 = (15/2)A3i? -7 and used the values k 2 — 0.17 
and fc 3 = 0.06 from Ref. [T7]. 

Nonlinear hydrodynamic corrections (~ 0.1%). 
The leading nonlinear hydrodynamic corrections 
are characterized by the coupling coefficient a/ A 3 
in the action. The size of this parameter can 
be estimated by comparing the Newtonian k 2 to 
the coupling constants in Lai's ellipsoidal models 
(e.g. Table 1 of [54]) to be uj 2 a/\ ~ 2 x 10~ 3 . 



The nonlinear self-coupling term in Eq. (A. 5) is 
smaller than the leading I = 2 term by a factor 
-285aw 2 /(572A) - 0.001. 

5. Spin corrections (< 0.3%). 

Fractional corrections to the tidal signal due to spin 
scale as 



5^ 



spin 



tidal 



"spin 



(A.7) 



where w max is the maximum rotational frequency 
the star can have before breakup, which for most 



NS models is > 2tt(1000 Hz). The observed NS- 
NS binaries which will merge within a Hubble time 
have spin periods of ~ 23 — 104 ms, and near 
the coalescence they will have slowed down due to 
e.g. magnetic braking, with final spin periods of 
~ 50 — 130ms. The fractional corrections to the 
tidal signal due to the spin are then < 0.3%. 

If the stars have spin, there will also be a spin- 
induced correction to the phase. As discussed 
above, the slow-rotation limit is likely to be the 
relevant regime for our purposes, and using sim- 
ilar methods as for the tidal corrections leads 
to a phase correction which scales as 8^ s cx 
n 2 R 2 /(r]M 2 x 1 / 2 )uj' 2 pin /(mi/R 3 ), where n 2 is the 
rotational Love number, which for Newtonian stars 
is the same as the tidal Love number k 2 and w S pm 
the spin frequency. The scaling of the spin term 
as tx x -1 / 2 shows that only at large separation do 
spin effects dominate over tidal effects, which scale 



as cx a; 



5/2 



Nonlinear response to the tidal field (< 3%). 

We have linearized in A. Including terms 



oc 
Eq. 



A 



gives 
(|A.5b of 



-4.8 x 10' 
-0.31fc 2 . 



li 



a fractional correction in 
-(83/7488)fc 2 i? 5 a; 5 / 2 /m 5 



fc 2 (m/M Q )- 10 / 3 (i?/km) 5 (//Hz) 5 / 3 



7. Viscous dissipation (negligible). 

There have been several analytical and numerical 
studies of the effect of viscosity during the early 
part of the inspiral, e.g. 0. They found that 
viscous dissipation is negligible during the early in- 
spiral if the volume-averaged shear viscosity ?7 s hear 
is 



Jfehear < l^' 



(- 



g cm 1 s . 



(A. 



The expected microscopic viscosity of NSs is 



10 



22 



10 14 g cnr 



9/4 / T V 
\10 6 KJ 



g cm *s , 



(A.9) 

which is orders of magnitude too small to lead to 
any significant effect. A variety of other likely 
sources of viscosity, e. g. the breaking or crum- 
pling of the crust, are also insignificant [S] [3] in 
the regime of interest to us. 

Thus, systematic errors in the measured value of A 
due to errors in the model should be O(10%), which is 
small compared to the current uncertainty of an order of 
magnitude in A. 
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